Electron-induced stabilization of ferromagnetism in Gai^Gd^N 
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Using ab initio band structure calculations and symmetry arguments, we show that the magnetic 
property of Gai-^Gd^N is drastically different from TM-doped GaN. The coupling between Gd 
atoms in the alloy is antiferromagnetic, but the ferromagnetic phase can be stabilized by introduc- 
ing electrons. Furthermore, we propose a model that may explain the recently observed colossal 
magnetic moments in this system, based on the polarization of donor electrons. 

PACS numbers: 75.50.Pp, 71.55.Eq, 71.70.-d, 71.20.Eh 
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INTRODUCTION 

The development of future spintronic devices will likely 
rely heavily on engineering materials that can have ef- 
ficient spin injection into semiconductors at high tem- 
perature P, U Among the most promising candidates 
for this task are some diluted ferromagnetic nitrides and 
oxides, which have small separation between mag- 
netic ions and small spin-orbit coupling at the band edge. 
Most of the materials that have been studied are semi- 
conductors doped with partially filled 3d transition met- 
als (TMs). Because of the strong coupling between the 
magnetic ions' 3d states and the coupling to the host p 
states, diluted magnetic nitride and oxide semiconduc- 
tors have been predicted and, in some cases, observed to 
show hole-mediated ferromagnetic behavior above room 
temperatures Q, 0, 0. However, high T c behavior of 
the 3d TM-doped nitrides and oxides is often impeded 
by the formation of precipitates Q or compensated by 
n-type carriers that are intrinsic for these wide gap semi- 
conductors 0, 0. Recently, nitrides and oxides doped 
with partially filled 4f rare-earth (RE) ions have been 
proposed as an interesting alternative to achieve high T c 
in these materials. Indeed, ferromagnetism has already 
been observed in Ga^u-jjGd^N, with Curie temperature 
larger than 400 K 0,|lQ. 

Despite the partial success, the nature of the host- 
impurity couplings, and of the ferromagnetic interactions 
in Gai-^Gd^N and other RE-doped nitrides and oxides, 
is not very well understood. On one hand, 4/ orbitals are 
more localized, thus the direct coupling between the 4/ 
ions is expected to be week. On the other hand, 4/ RE 
elements can have larger magnetic moments than the 3d 
elements, and, unlike the d states, / electrons can couple 
strongly with the host s electrons, leading to the possi- 
bility of electron-mediated ferromagnetism in these ma- 
terials. Furthermore, colossal magnetic moments of Gd 
in GaN have been observed, which has a close connection 
to the observed ferromagnetism in this system But 
the origin of the observed colossal magnetic moments is 
not known. 

To get a deep understanding of the nature of mag- 
netic coupling in these RE-doped systems, in this work, 



we perform total energy and band structure calcula- 
tions of Gai-^Gd^N for a few concentrations, including 
the hypothetical GdN in the zinc-blende structure. We 
show that the interactions between Gd atoms in undoped 
Gai-^Gd^N are antiferromagnetic in nature, but this 
can be changed by introducing donors to the material; 
different from most 3d diluted magnetic semiconductors, 
ferromagnetism in Gai-^Gd^N is electron-mediated [loj . 
because the splitting of the conduction band, induced by 
the s-f coupling, is larger. 



METHODOLOGY 

The total energy and band structure calculations in 
this study were performed using the general poten- 
tial linearized augmented plane wave (LAPW) method 
(l"l| within the density functional theory and the local 
spin density approximation (LSDA) for the exchange- 
correlation potential. The muffin-tin (MT) radii are 1.48 
Bohr for N, 2.50 Bohr for Gd, and 2.05 Bohr for Ga. The 
Gd bp and Ga 3d semicore states are treated in the same 
footing as the other valence orbitals. The Brillouin-zone 
integrations were performed using a 4x4x4 Monkhorst- 
Pack special /c-points [l2T | and equivalent Appoints for 
the superstructures. We find that the LSDA-calculated 
lattice parameters for both GaN and GdN are in excel- 
lent agreement with experiment. Our calculated lattice 
parameters for GaN and GdN in the zinc-blende struc- 
ture are 4.497 and 5.291 A. The difference in the lattice 
constant between ferromagnetic (FM) and antiferromag- 
netic (AFM) configurations of GdN is less than 0.2%. It 
should be noted that pure GdN usually crystallizes in 
the NaCl structure. We calculated the lattice parameter 
of GdN in this structure, finding it to be 4.978 A, which 
is in excellent agreement with the experimental value of 
4.974 A0|. For the alloys, we assume that Vegard's rule 
is obeyed for the lattice constant, and the internal atomic 
positions are fully relaxed by minimizing the total energy 
and quantum mechanical forces. In this work we study 
Gd incorporated in zinc-blende GaN and assume that the 
physical properties are similar for Gd incorporated in the 
more stable wurtzite GaN. 
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FIG. 1: (Color online) Band structure of the ferromagnetic 
GdN in the hypothetical zinc-blende structure. The labels 
indicate the majority character of the band at the T point. 




FIG. 2: Wavefunction square of the (a) spin-up VBM, (b) 
spin-down VBM, (c) spin-up CBM, and (d) spin-down CBM 
at the zone center. 



RESULTS AND DISCUSSION 

To understand the effect of Gd incorporation in GaN, 
we first look at the band structure (Fig. 1) of the hypo- 
thetical zinc-blende ferromagnetic GdN compound and 
compare it with that of MnN Q. We find that in the 
tetrahedral environment, the crystal field splits the / or- 
bitals into t\, t%, and aq states, which can couple to other 
s, p, and d states with the same crystal symmetry. For 
Gd, the majority spin-up 4/ states are fully occupied and 
are below the valence band maximum (VBM) state at T 
with mostly N 2p character and £2 symmetry. The mi- 
nority spin-down 4/ states are fully unoccupied and are 
above the conduction band minimum (CBM) at T. Un- 
like the 3d TM, the Gd 5d bands, with e and ti symmetry, 
are high in energy, and thus arc fully unoccupied in the 
conduction band. The unique position and characters of 
the Gd 4/ and 5d orbitals lead to some unusual features 
of GdN that are very different from those of Mn-doped 
GaN. (i) The VBM state at T, with t2 P character, couples 
strongly with the Gd tid and £2/ states. In the spin-up 
channel, this state is pushed up by the Gd ti / state below 
and pushed down by the Gd tid state above. In the spin- 
down channel, this state is pushed down by both the Gd 
t2f and tid states. However, the p-d coupling is larger 
in the spin-up channel, because the potential exchange 
lowers the spin-up Gd tid orbital energy; the net effect is 
that the exchange splitting at the VBM (Nq/3 ~ -0.01 
eV) is much smaller than that observed in TM-doped 



GaN. Furthermore, because of the strong p-d coupling at 
F, which pushes the ti v state down, the top of the va- 
lence band is not at T as for MnN 0. (ii) Unlike 111 
TM-doped GaN, where s-d coupling is not allowed in Td 
symmetry, s-f coupling is allowed. Because the CBM at 
r has the s character, in the spin-up channel, it is pushed 
up by the occupied Gd a\f state, and in the spin-down 
channel, it is pushed down by the unoccupied Gd a\f 
state. This large kinetic exchange leads to a large nega- 
tive spin exchange splitting near the CBM (iVoa ~ —0.4 
eV). In contrast, the potential exchange induced splitting 
in TM-doped GaN is smaller and positive. The wave- 
function square of the VBM and CBM states at T are 
plotted in Fig. 2. We see that the spin-up VBM shows a 
large antibonding p-f character. The / character in the 
spin-down VBM is relatively small because the Gd / and 
VBM energy difference is large in the spin-down chan- 
nel. For the CBM states, we find that the spin-up CBM 
shows an antibonding f-s character because the spin-up 
/ state is below the CBM, whereas the spin-down CBM 
shows a larger bounding f-s character because the unoc- 
cupied spin-down / level is above and closer to the CBM. 
These results are consistent with the energy position and 
symmetry arguments above. 

For diluted GaN:Gd magnetic semiconductors, the 
general features are the same as in pure GdN. The ma- 
jority Gd 4/ orbitals of Gd are always inserted below 
the VBM of GaN, whereas the minority Gd 4/ orbitals 
are always above the CBM. The Gd 5d orbitals, on the 
other hand, always have higher energy, above the CBM 
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of GaN. These features can be seen in Figs. 3 and 4. In 
Fig. 3 we show the band structure for Gai-^Gd^N with 
x = 0.25. Figure 4 shows the density of states (DOS) 
for Gai-^Gd^N with x — 0.0625. We observe that the 
spin-up Gd 4/ orbitals are resonant within the valence 
band with large band widths, whereas the spin-down Gd 
4/ orbitals are more localized. We can also observe that 
the DOS at the VBM is very high, indicating a very large 
valence effective mass in this material. The large valence 
effective mass comes from the coupling between the VBM 
and the Gd d orbitals above, which can be seen in Fig. 
4. The repulsion between these orbitals pushes the VBM 
at r down, increasing the effective mass. This indicates 
that this material will have a very small hole mobility. On 
the other hand, the DOS of the conduction band is very 
small. This can be understood by the repulsion between 
the 4/ levels and the CBM. This interaction pushes the 
CBM at T down in the spin-down channel, decreasing the 
electron effective mass. We find that, again, both NqOc 
and NqP are negative for Gai-^Gd^N, with the magni- 
tude of N a much larger than that of N (3, as can be 
seen in Fig. 3. This effect is an indication that, for 
Gai_ x Gd x N, the stabilization of ferromagnetism will be 
more efficient when the carriers are electrons instead of 
holes. 

Because of the highly localized character of the / or- 
bitals, a natural question is whether LSDA can suc- 
cessfully describe the system. Fortunately Gd has a 
4/ 7 5d 1 6s 2 valence configuration, thus it is isovalent with 
Ga. The majority-spin 4/ orbitals are completely filled, 
whereas the minority-spin 4/ orbitals are completely 
empty. This situation is similar to Mn substituted II- VI 
diluted semiconductors . The fact that the calculated 
lattice parameter for GdN is in good agreement with ex- 
periment indicates that our LSDA description is reason- 
able. To further test the adequacy of the LSDA, we have 
performed general gradient approximation (GGA) [lfij ] 
and model LDA+U type calculations flij. We find that 
besides small quantitative changes (Table I), the qualita- 
tive picture described in this paper is unchanged. This 
can be understood because GGA and LDA+U does not 
change the /-electron occupation in this system. 

Our band structure calculation show that undoped 
zinc-blende GdN or Gai-^Gd^N alloys ml, are semi- 
conductors because Gd is isovalent to Ga 18] . Our total 
energy calculations find that undoped GdN and GaGdN 
in the zinc-blende phase are more stable in the antifer- 
romagnetic phase. The calculated results are shown in 
Table I for several Gd concentrations x and configura- 
tions. There are always two Gd atoms in the unit cell. 
These results are expected because no charge carriers are 
created in this system, therefore, the spin-splitting at the 
band edge present in the FM phase will not lead to any 
energy gain compared to the unsplit band edge of the 
AFM phase. On the other hand, the superexchange 
coupling between the Gd / states, mediated by the N p 
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FIG. 3: (Color online) Band structure for the ferromagnetic 
configuration of Gai-^Gd^N with x — 0.25. 



orbitals, stabilizes the AFM phase. However, compared 
to TM-doped GaN, the energy differences (AEfm-afm) 
between the FM and AFM states for these systems are 
small, which is due to the highly localized character of 
the / orbitals. 

In order to stabilize the FM phase, it is necessary to 
insert carriers into the system. As GaN is usually n-type 
doped with free electrons in the conduction band 0, we 
have simulated the effect of donors by introducing elec- 
trons into Gai-^Gd^N. In the AFM configuration, the 
CBM of each spin channel have the same energy, there- 
fore, they will be occupied in both spin channels. In the 
FM configuration, because of the strong s-f coupling, the 
CBM has negative spin-exchange splitting, i.e., the spin- 
down CBM has a lower energy. Consequently, when elec- 
trons are added, they will prefer to go to the spin-down 
CBM of the FM phase, stabilizing the FM configuration. 
This expectation is supported by our direct total energy 
calculations of charged systems, where the FM phase be- 
comes more stable. We have considered two different 
configurations for x = 0.125, where the Gd atoms are 
nearest neighbors and next-nearest neighbors. We find 
that the effective exchange interaction between the two 
Gd [Jjj(R)] decreases as the distance between impurities 
becomes larger. 

Several experimental studies 0, 11 @ of Gai-zGd^N 
claim that ferromagnetism is observed up to room tern- 
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FIG. 4: (Color online) Density of states of Gai-^Gd^N with 
x — 0.0625. The blue (dashed) curves are partial DOS of Gd 
/ levels and the red (dotted) curves are partial DOS of Gdd. 
For clarification, the intensity of the DOS related to the Gd 
d levels is multiplied by a factor of three. 



TABLE I: Energy differences (AEfm-afm) between 
ferromagnetic and antiferromagnetic configurations in 
Gai-zGd^N. Energies are in meV; nn and nnn refer to 
the two Gd atoms as nearest and next-nearest neighbors, 
respectively, in the unit cell. 



x Charge (electrons/Gd) AEfm-afm 

1 25.0 

1 22.0 (GGA) 

0.25 2.6 

0.125 (nn) 4.4 

0.125 (nn) 0.1 -28.4 

0.125 (nnn) 1.3 

0.125 (nnn) 0.1 -19.5 



perature. Our calculations show that these observed FM 
phases should be related to the n-type character of the 
Gai-^GdajN samples. This is supported by photolumi- 
nescence measurement, which shows the presence of large 
defect concentration in the sample with FM ordering [20j . 
A recent measurement also shows that the presence of 
FM is correlated with Oxygen, a known donor defect in 
the sample with a concentration of about 10 18 /cto 3 . 
These observations are consistent with our model. 

One of the unique features observed in Gd-doped GaN 
is the colossal magnetic moment Q. It was observed 
that, at very low Gd concentrations (~ 10 16 /cm 3 ), the 
total effective magnetic moment per Gd atom could be 
as high as 4000 /j,b- This is quite unusual because the 
atomic moment of Gd is only 8 [Lb , and in the hypothet- 
ical zinc-blend GdN the moment per Gd is only 7 hb- 
When the Gd concentration increases, the effective mo- 
ments per Gd decrease. This observation was explained 



through a phenomenological model, which assumes that 
Gd atoms polarize the GaN matrix in a certain radius 
around them and each host atom within the radius has 
an induced magnetic moment pq. Although this model is 
very interesting and fits experimental results quite well, 
at least in the low-concentration region, it does not pro- 
vide a mechanism of the induced magnetic moment. Our 
band structure calculations above, however, provide a 
possible explanation of the observed huge magnetic mo- 
ments: When Gd is introduced into GaN, each Gd in- 
troduces unoccupied / states above the spin-down CBM. 
The symmetry-allowed coupling between the / states and 
states below it can create localized levels with energy be- 
low the original CBM. When the system also contains 
enough shallow donor impurities, as reported in the ex- 
periment, two things will occur: (i) the donor electrons 
of the impurities will occupy the Gd-induced empty lev- 
els, if the donor levels of the impurities are higher than 
the Gd-induced level; (ii) consequently, the system will 
become ferromagnetic. When enough levels are created 
by Gd below the donor level, which is possibly due to 
the symmetry-allowed s-f coupling, the induced mag- 
netic moment per Gd atom can be very high in the ex- 
treme impurity limit when the Gd-Gd distances are very 
large. Translating this into a real-space representation, 
it is equivalent to say that atoms with a certain radius 
ro around each Gd atom has its spin-down energy levels 
below the donor levels, and thus they are polarized by 
the Gd atom with opposite spins. As Gd concentration 
increases, the spheres surrounding each Gd atom start to 
overlap. Although some more levels could be pushed be- 
low the donor level, thus increasing slightly the effective 
radius ro , eventually the overlap will reduce the effective 
volume affected by each individual Gd, and therefore, re- 
duce the magnetic moment per Gd, as observed in exper- 
iment. However, unlike the previously proposed model, 
our model has two distinct physical characteristics: (i) To 
have a large induced magnetic moment, the system must 
have enough donor impurities; (ii) Our model suggests 
that the induced magnetic moment has the opposite sign 
of the Gd magnetic moment. This indicates that when 
the Gd concentration increases, or the donor carrier den- 
sity decreases, the magnetic moments in the system will 
change sign. Experimental testing of our model is called 
for. 

It is interesting to point out that this electron-induced 
stabilization of ferromagnetic order and large magnetic 
moments may also occur in TM-doped semiconductors 
flfT |. This is because, despite that s-d coupling is not 
allowed under Td symmetry, it will be allowed when the 
symmetry of the system is reduced either by strain or by 
chemical ordering. For example, under (001) strain, the 
ai(s) state can couple to the split ai(e,j) state, whereas 
under (111) strain or in wurtzite crystal, the ai(s) state 
can couple to the split a\(td) and a\{t p ) state. However, 
in general, the s-d coupling is expected to be smaller than 



5 



the symmetry-allowed s-f coupling in these tetrahedral 
semiconductor systems. 



SUMMARY 

In summary, we have investigated in detail the elec- 
tronic and magnetic properties of Gai-^Gd^N. We find 
that because of the coupling between the Gd / and host 
s states, this system shows some unique behavior that is 
drastically different from TM-doped GaN. The exchange 
splitting at the CBM is negative. The coupling between 
Gd atoms is found to be antiferromagnetic, if no donors 
are present, but it become ferromagnetic when enough 
donors are present in the system. We also proposed a 
model that may explain the colossal magnetic moments 
observed in this system in the very dilute limit, showing 
that it should be directly related to the polarization of 
donor electrons. 
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